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We use density functional theory and the van der Waals density functional (vdW-DF) method to 
determine the binding separation in bilayer and bulk graphane and study the changes in electronic 
band structure that arise with the multilayer formation. The calculated binding separation (distance 
between center-of-mass planes) and binding energy are 4.5-5.0 A (4.5-4.8 A) and 75-102 meV/cell 
(93 — 127 meV/cell) in the bilayer (bulk), depending on the choice of vdW-DF version. We obtain the 
corresponding band diagrams using calculations in the ordinary generalized gradient approximation 
for the geometries specified by our vdW-DF results, so probing the indirect effect of vdW forces on 
electron behavior. We find significant band-gap modifications by up to -1.2 eV (+4.0 eV) in various 
regions of the Brillouin zone, produced by the bilayer (bulk) formation. 
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I. INTRODUCTION 

Selective modification of band gaps (band-gap engi- 
neering) by atomic-scale design of materials is a pow- 
erful concept in electronic and photonic development. 1 
Band gaps can be altered by, for example, introducing 
dopants, defects or by exploiting finite size effectsP^The 
physical origin of band-gap variations is a modification 
in the charge distribution in concert with wavefunction 
hybridization and modification. 

Dispersive or van der Waals (vdW) interactional also 
alter the distribution of electronic charges and hence the 
electron band- structure. A direct effect is evident, for 
example, by considering the formation of the double- 
dipole configuration 5 which is the electrostatic signature 
and inherent nature of a pure vdW binding P*^ In addi- 
tion, there are also indirect, geometry -induced effects of 
vdW binding on electron behavior. These indirect ef- 
fects arise when two material fragments come close to 
one another, thereby changing the local electron environ- 
ments as compared to isolated fragments. For example, 
vdW-binding can cause smaller amounts of net charge 
transfer within individual vdW-bonded fragments. 6 Also, 
wavefunction hybridization will certainly arise when ma- 
terial fragments approach one another, even if this hy- 
bridization does not significantly contribute to the bind- 
ing itself (in purely dispersive interaction). Wavefunc- 
tion hybridization and Pauli exclusion cause scattering 
of the surface-state electrons in physiosorption of acenes 
and quinones on Cu(lll), even if there is no net charge 
transfer 

The strength of vdW interactions in surface/adatom 
systems (including adclusters or complete overlayers) 
can be controlled by the precise choice of the sur- 
face materiaP^ and by the surface morphology (flat, 
stepped, pyramidal, etc.). 12 It is therefore important 
to quantify the extent to which van der Waals (vdW) 
bonding can modify the electron dispersion, that is, the 
band structure. Such a study is now possible, since re- 
cent development of the vdW density functional (vdW- 



DF) methocP^n^ enables systematic (theoretical) ex- 
plorations of bonding in sparse materials within density 
functional theory (DFT)P 

In this paper we investigate indirect, geometry-induced 
effects of vdW bonding on electronic structure. In partic- 
ular we separately study the effects of (1) the enhanced 
charge density in regions where the tails of vdW-bonded 
material fragments overlap, see left panel of Fig. [I] (2) the 
hybridization of wave functions, and (3) the redistribu- 
tion of charge density due to mechanisms that are not 
inherent to the vdW interaction (such as local displace- 
ments due to electrostatics and Pauli-repulsion) , see right 




FIG. 1: (Color) Charge density in vdW-bonded graphane 
bilayers. The left panel shows a schematics of the non- 
selfconsistent bilayer (BL) charge density obtained by adding 
two monolayer (ML) densities and resulting into a charge 
enhancement in the inter layer region. The right panel de- 
picts calculated self-consistent (on the GGA level) charge re- 
arrangements. Our Bader analysi s 23 * 24 * shows that the charge 
rearrangements are not accompanied by any significant net 
charge transfer into the middle region. Color coding: light 
large spheres represent C atoms, black small spheres repre- 
sent H atoms; accumulation of charge is indicated by the pink 
(light) isosurface, depletion by the purple (dark) isosurface. 



panel of Fig. [T] 

We focus on the band stru cture of layered systems of 
the macromolecule graphane j 16 * 17 * a fully hydrogenated 
derivative of graphene.^ The top panel of Fig. [2] shows 
the atomic structure of (the stable chair conformation 
of) monolayer (ML) graphane, consisting of a (slightly 
buckled) graphene backbone with H atoms attached in 
alternating fashion above and below the carbon plane. 
Bilayer graphane with possible high-symmetry structure 
shown in the bottom panels of Fig. [2] present a possible 
new system i n addition, we also include a study of 
a possible bulk-graphane crystal. 

Graphane adds to the wealth of carbon-based materi- 
als that are considered as pr omisin g materials for near- 
future nanoelectronic devices ! 21 * 22 * Pure graphene has a 
zero band gap and extraordinary conduction properties. 
Electronic devices, however, also require semiconducting 
and insulating materials. Such materials can be obtained 
from pure graphene as derivatives either in the form of 
graphene nanoribbons 25 26 or by chemical modification 
through adsorbates! 27 | 2S | Monolayer graphane belongs to 
the last-mentioned group of derivatives. DFT calcula- 
tions predict a large band-gap semiconductor naturej^ 
the more advanced GW methocP^ predicts an insulat- 
ing naturepS Graphane has been proposed theoretically 
to serve as a natural host for graphene quantum dots^ 
or graphene nanoribbons for nanoroads. 32 Furthermore, 
doped graphane has been recently predicted to be a high- 
ly superconductor P^l Such potential application of the 
graphane structure makes it interesting to explore pos- 
sibilities to (locally) modify the the electronic behavior 
either by selective hydrogen removaP^^or by geometry- 
induced band-structure modifications. 

The paper is organized as follows. In Sec.|TI| we give a 
survey of all considered high-symmetry graphane bilayer 



configurations. Section [HI] presents our computational 
method. In Sec. IV we present, analyze and discuss our 
results. Section I VIII summarizes our work and contains 
our conclusions. 



II. HIGH-SYMMETRY GRAPHANE BILAYERS 

The set of lower panels in Fig. [2] shows all six high- 
symmetry arrangements of bilayer (BL) graphane. These 
can be grouped into two different types. In a- type BL, 
the graphane sheets are interlocked with each other. In 
/3-type BL, the H atoms from different graphane sheets 
(located between the sheets) sit on top of each other. We 
calculate and compare all of these configurations that 
make up the a- and /3-type sets of stacking configura- 
tions. 

The set of different (high-symmetry) arrangements for 
the BL systems are found as follows. We label the sheets 
according to the location of vacancy in the C backbone 
in the unit cell (A, £?, or C sites). In addition, the dis- 
tortion of the C backbone along the z-direction (+ or 
— ) of the first occupied C site (counted along the main 
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FIG. 2: Structure of monolayer (top panel) graphane and of 
high- symmetry graphane bilayer (BL) configurations. The la- 
beling (A+A_, . . . ) provides an unambiguous and exhaustive 
identification of all high-symmetry atomic configurations for 
the BL system as explained in the text. Light large spheres 
represent C atoms, black small spheres represent H atoms. 



diagonal starting from the vacancy stacking) is indicated 
as a subscript label. In all BL, the first layer can be ar- 
bitrarily chosen to be an A+ layer. The second layer is 
placed on top of the first one of the following actions: 
(i) copying the bottom layer and moving it along the 
^-direction (A + A + ); (ii) flipping the bottom layer and 
moving it along the z-direction (A + A_); (hi) as in (i) 
and additionally moving it along the long diagonal of the 
2D graphane lattice by one third (A + £? + ); (iv) as (ii) 
and additionally moving it along the long diagonal of the 
2D graphane lattice by one third (A + £?_); (v) as in (i) 
and additionally moving it along the long diagonal of the 
2D graphane lattice by two thirds (A + C + ); (vi) as in (ii) 
and additionally moving it along the long diagonal of the 
2D graphane lattice by two thirds (A + C_). 



III. COMPUTATIONAL METHOD 

A. vdW binding 

We map out the energy variation of bilayer (bulk) 
graphane as a function of the separation between two 
graphane sheets (the c-parameter of the bulk unit cell) 
using (non-selfconsistent) vdW-DF calculations. For 
BL, we employ supercells with our optimized lxl 
graphane in-plane lattice parameters [a\ = (a, 0,0), <22 = 
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(a/2, >/3/2 • a,0) with a = 2.532] and a height of 30 A. 
For the bulk, we optimize the c-parameter of the periodic 
unit cell [as = (0, 0, c)] starting from the optimal value 
of the BL separation. 

Our calculations combine selfconsistent DFT calcula- 
tions in the generalized gradient approximation (GGA) 
with three (non-selfconsistent) versions of the vdW-DF 
method. The GGA calculations are performed with the 
planewave pseudopentiaP^ code Dacapoj^ using PBEP^ 
for exchange and correlation. We use a planewave cut- 
off of 500 eV and a 4x4x1 (4x4x2) k-point sampling^. 
The three versions of the vdW-DF method that we use 
are (i) the nonlocal correlation functional of Dion et 
a/P^ in conjunction with revPBEP^ for exchange (vdW- 
DF1), (ii) the same correlation functional but with the 
exchange part of the C09 functional 15 (vdW-DFl-C09 x ), 
and (hi) the most recent version of the vdW-DF method, 
Ref. [Ml (vdW-DF2). The latter version uses the refitted 
form of the PW86 functional (rPW86 x ) 39 for exchange. 
We obtain total energies as 
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Here, E^ l [n] is the energy obtained from one of the non- 
local functionals of Refs. [T3l and [T4| and Eo[n] is given 
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where i^ywN is the VWN-LDA 40 correlation energy and 
the subscript V denotes the version of the exchange func- 
tional (revPBE x , C09 x , or rPW86 x ). We define the layer 
binding energy as 



^bind(^cmp) = ^vdW-DF(4mp) ~ ^vdW-DF(^c 



OO 



(3) 



Here, d cmp is the distance between the center-of-mass 
planes in each graphane sheet of the monolayer. 

Our numerical evaluation of Eq. ([3| proceeds in the 
same way as described in Refs. [71 and |41 44 . In par- 
ticular, because of a small but nonnegligible sensitivity 
of the nonlocal correlation on the exact positioning of 
atoms with respect to the density grid, we avoid a di- 
rect comparsion of E^ l [n] for BL configurations with dif- 
ferent ML separations. Instead we evaluate the layer- 
binding energy by comparing changes in the nonlocal 
correlation arising between the actual configuration and 
a reference that keeps the same alignment of atoms and 
grid points. Specifically, for all configurations we eval- 
uate the change in nonlocal correlation as AE^ l [n] = 
pq N - Ef, p N " £ c nl Q N • Here E$ PQ [n] is the non- 
local correlation energy of the full BL configuration (with 
one ML in P and one in Q) and E^ p [n] (E^ 1 q[ti\) is the 
nonlocal energy of the configuration where one ML has 
been removed from Q (P) while the other is kept at pre- 
cisely the same location P (Q) as in the BL configuration. 
Further details on our approach to increase the accuracy 
of vdW-DF are provided in Refs. ITlandHTI 
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FIG. 3: Calculated layer-binding energy variations as func- 
tions of the center-of-mass (cmp) separation d cm p between two 
graphane sheets. The top panel compares energy variations 
of the six different high-symmertry graphane configurations 
calculated with vdW-DFl-C09 x . The a- type bilayer have in 
general lower energy than /3-type configurations. The A+A+ 
on-top configuration shows the strongest binding. The bot- 
tom panel compares the different version of vdW-DF for the 
A+A+ configuration. All versions of vdW-DF predict quali- 
tatively the same energy variations but the detailed numerical 
values of the binding separations and energies vary. The in- 
sert emphasizes that GGA calculations provide no meaningful 
account of the binding. 



B. Band structure 

We determine band structures with pure GGA calcu- 
lations for various k points. We fix the BL separation 
(bulk c-parameter) to the value calculated with vdW- 
DF, so probing indirect effects of vdW binding on elec- 
tron behavior. The (selfconsistent GGA) input density 
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cfcmp [A] 

E hind [meV/cell] 
£ gap at T [eV] 
£ gap at K [eV] 
E gap at M [eV] 



bilayer graphane 
vdW-DFl vdW-DF2 vdW-DFl-C09 x 



5.0 

78 

3.46 (-0.08) 
11.86 (-0.30) 
10.40 (-0.46) 



4.75 
77 

3.53 (-0.01) 
11.75 (-0.41) 
10.31 (-0.55) 



4.5 
101 

3.61 (+0.07) 
11.61 (-0.55) 
10.24 (-0.62) 



bulk graphane 



c[A] 

^bind 



-oina [meV/cell] 
E gap at T [eV] 
E gap at K [eV] 
E gap at M [eV] 
E gap at A [eV] 
E gap at H [eV] 
E gap at L [eV] 



vdW-DFl 
4.8 
93 

7.23 (+3.69) 
12.77 (+0.61) 

10.00 (-0.86) 
3.52 (-0.02) 

11.63 (-0.53) 
11.27 (+0.41) 



vdW-DF2 
4.7 
94 

7.50 (+3.96) 
12.83 (+0.67) 
9.95 (-0.91) 
3.62 (+0.08) 
11.58 (-0.58) 
11.27 (+0.41) 



vdW-DFl-C09 x 
4.5 
127 

8.01 (+4.47) 
12.66 (+0.50) 
9.85 (-1.01) 
3.88 (+0.34) 
11.48 (-0.68) 
11.28 (+0.42) 



TABLE I: Binding separations d cmp (the c- lattice param- 
eter) layer-binding energies energies i^bind, and band gaps 
E gap at several k points in the Brillouin zone of bilayer (bulk) 
graphane. All quantities are calculated with three versions of 
vdW-DF. For the band gaps, the difference with respect to a 
graphane monolayer are given in parentheses [a negative value 
corresponds to a decreased band gap in the bilayer (bulk)]. 



for these (non-selfconsistent GGA) band-structure calcu- 
lations is obtained using a planewave cutoff of 500 eV and 
a 20x20x1 (20x20x5) k-point sampling 37 for the bilayer 
(bulk). 

For the bilayer, the Brillouin zone (BZ) is two- 
dimensional and relevant k points are T — (0,0,0), 
K = (2/3,1/3,0) and M = (1/2, 1/2, 0)P All special 
points are given in units of the reciprocal lattice vectors. 
We calculate the band variations along the paths KT, 
rM, and KM. 

For the bulk, the BZ is three-dimensional. Therefore 
non-zero values of k z are important and we also calculate 
the band variations along HA, AL, and LH. Here, the 
special points are A = (0,0,1/2), H = (2/3,1/3,1/2), 
and L = (1/2,1/2,1/2)™ 
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FIG. 4: Electronic band structure of bilayer (BL) and bulk 
graphane. The top panel shows the overall band diagram 
along KT, FM, and MK of a A+A+ BL at the binding sepa- 
ration predicted by vdW-DF2. Apart from the fact that each 
band exists as a pair of bands, the BL band-structure quali- 
tatively agrees with that of of a monolayer (ML, see Fig. 3 in 
Ref. I16p . The bottom panel shows the /c-dependent (direct) 
band-gap variation in the BL and bulk with respect to the 
band-gap variation in the ML. For the bulk, we also show the 
band-gap variation along HA, AL, and LH. In the BL, at and 
around the F point, the gap is essentially unchanged. How- 
ever, away from the T point significant band-gap reduction 
is observed. In the bulk, local band-gap modifications can 
be both positive and negative and their absolute can be even 
larger than in the BL. 



IV. RESULTS: PREDICTED PROPERTIES OF 
BILAYER AND BULK GRAPHANE 

Figure [3] shows the calculated variations in layer- 
binding energies as functions of the separation between 
the center-of-mass planes ('cmp') of the two graphane 
monolayers (ML) in a bilayer (BL). The top panel com- 
pares the energy variations for the configurations with 
different stackings using vdW-DFl-C09 x . The energy 
variations split according to the grouping into a- and 
/3-type configurations. The a- type configurations have a 
smaller binding separation and a higher binding energy; 
the stacking shows the strongest bonding.™ 

The bottom panel compares the energy variations of 
the A + A + BL for vdW-DFl, vdW-DFl-C09 x , and vdW- 



DF2. Qualitatively, all functionals yield the same energy 
variations. The insert shows the energy variation for the 
A+A+ configuration obtained from pure PBE calcula- 
tions and illustrates that no meaningful binding is pre- 
dicted without an account of vdW forces. 

Table [T] lists and compares numerical results for the 
calculated binding separations and layer-binding energies 
for the BL. The binding separations and energies range 
from 4.5 A to 5.0 A and 75 meV/cell to and 102 meV/cell, 
depending on the version of vdW-DF.™ The binding en- 
ergy is comparable to that in a graphene BLf^ (94 meV 
using vdW-DFl). 

Table U also lists the calculated lattice constant c and 
the corresponding layer-binding energies for a fictitious 
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bulk crystal of graphane. It is possible that such a 3D 
graphane system might eventually be synthesized. We 
here present predictions of the expected structure, using 
our analysis of stacking in the BL as starting point. In 
particular we focus on A+A+ stacking and find that the 
lattice constant essentially coincides with d cmp in the bi- 
layer. The binding energy is slightly increased in the bulk 
and varies between 93 eV and 127 meV. These numbers 
also compare to the binding in graphite^ (100 eV using 
vdW-DFl). 

The top panel of Fig. [4] presents the overall PBE band- 
diagram for the stacked graphane BL at the bind- 
ing separation predicted by vdW-DF2. Corresponding 
band diagrams at vdW-DFl or vdW-DFl-CO x binding 
separations are qualitatively similar. Apart from the fact 
that each band occurs as a pair of bands, the band struc- 
ture also agrees qualitatively with that of the ML (see, 
for example, Ref . [T6]) . 

The bottom panel of Fig. [4] summarizes some differ- 
ences between the BL (bulk) and the ML, document- 
ing the changes occurring in the ^-dependent band gap 
with the BL (bulk) formation. We plot the differ- 

ences Ai^ ulk (fc) = £ g B a L P /bulk (fc) - Ef^{k) along KT, 
TM, and MK. In addition, for the bulk, we also plot 
A£ W ( k ) = £gap k ( fe/ ) - ^ g M apW- Here, k' is along 
HA, AL, and LH in the bulk; k is the corresponding 
k point along KT, TM, and MK in the ML (and there- 
fore k x = k' x1 k y = k' y1 while k z = and k' z = 1/2). A 
summary of the numerical values of band gaps at the spe- 
cial points (calculated with various choices of vdW-DF 
and corresponding BL binding separations or bulk lattice 
constants c) and their deviations from the corresponding 
values in the ML is given in Table [l| 

We find large modifications of the band structure, indi- 
rectly induced by the vdW interactions and summarized 
by the /c-dependent band-gap differences. In the BL, the 
direct band gap can deviate by up to ~ 0.8 eV (between 
K and T) with respect to the ML gap (see bottom panel of 
Fig.|4|. In the bulk, deviations can be as large ~ — 1.2 eV 
(~ +4 eV) in some regions of the Brillouin zone (BZ) near 
the H point (r point). 



V. DISCUSSIONS: BILAYER GRAPHANE 

Focusing on the direct band gap in BL graphane, we 
find modifications that are strongly /^-dependent. At the 
K and M points (and in other regions), the modifications 
are significant. At the T point, where the gap is small- 
est in the ML (and in the BL), no modifications occur, 
rendering the BL system electronically similar to the ML 
system. Nevertheless, qualitative understanding of the 
origin of the different modifications in the various regions 
is important to gain further insight into the relevance of 
vdW interactions for materials band-structure. 

In the following, we explore the band-gap modifications 
upon formation of graphane BLs in more detail. We focus 
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FIG. 5: Band-origin of gap modifications, self-consistent 
charge rearrangements and wave function hybridization. The 
top and bottom panels show the band structure in the LCB 
and UVB respectively. ML graphane bands have been in- 
cluded to highlight that, in the BL, bonding and antibonding 
hybrid bands are formed. For the MLs, the band index n is 
annotated. Also, we have marked some combinations of bands 
and k points for which wavefunction will be analyzed later. 
The mid-panel shows the variation of direct band-gap differ- 
ences between the BL and ML (solid black), the variation of 
UVB energy differences (dashed-dotted black), and the varia- 
tion of LCB energy differences (dashed black) . We also show 
the variation of band-gap differences calculated from the non- 
selfconsistent BL charge density obtained by adding two ML 
densities (dashed light gray). 

on the indirect, geometry -induced effects of vdW binding 
outlined in Fig. [I] the effect of a modified electronic envi- 
ronment in the region between the MLs that form the BL 
arising from a superposition of two ML electron densities 
(see left panel); the effect of self-consistent charge rear- 
rangements (on the GGA level) of this superpositioned 
density (see right panel) ; the effect of potential hybridiza- 
tion of wave functions (WFs) (not shown in the figure). 

Our analysis suggests that the band-gap modifica- 
tions should be interpreted as a concerted interplay be- 
tween WF hybridization and electrostatic interaction be- 
tween the hybridized WFs with the modified environ- 
ment. The relevant WFs are unoccupied conduction- 
band (CB) WFs. Self-consistent (SC) charge rearrange- 
ments (with respect to the superposition of ML densities) 
do not play a significant role. 

A. Band-origin of gap reduction 

In the mid-panel of Fig. |5j we show that the main con- 
tribution to the observed band-gap variation comes from 
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the lower conduction-band (LCB). We plot the varia- 
tions of the upper valence-band (UVB) energies (dashed- 
dotted curve) and of the LCB energies (dashed curve), 



AE X xx(fc) = (£fxxW " £f L ) " (E^x(k) ~ E> 



(4) 



Here E^^ AL (k) * s tne ener gy of the highest (occupied) 
valence band (XXX = UVB) or the energy of the lowest 
(unoccupied) conduction band (XXX=LCB) at k in the 
BL/ML and E® L/ML is the Fermi level in the BL/ML. 
(The Fermi level is here defined by the energy of the 
highest occupied state). 

We find that |A^lcb(^)| is typically much larger than 
|A£'uvb(^)|- The latter is essentially zero. In the UVB, 
nonzero contributions to the band-gap modifications are 
only found in small regions. As can be seen from the bot- 
tom panel of Fig. [5] where we contrast BL VBs with ML 
VBs, the topmost BL UVBs are non-degenerate there, 
indicating WF hybridization. Nevertheless, in regions 
where \AE gap (k)\ is large, only |A£^lcb(^)| contributes. 
Thus, we assign the modifications of the band-gap varia- 
tion in the BL primarily to the modifications in the LCB 
energy variation. 



B. Self-consistent charge rearrangements 

In the mid-panel of Fig. [5] we also contrast band- 
structure modifications obtained from the SC BL charge 
density with band-structure modifications obtained from 
the non-SC charge density, constructed as a superposi- 
tion of ML densities. The solid black line (gap) corre- 
sponds to the SC case, the dashed light-gray line (gap*) 
corresponds to the non-SC case. At the displayed reso- 
lution, the curves cannot be distinguished. We find that 
the differences between both band-gap variations are at 
the meV level. Thus, the charge rearrangements^ shown 
in the right panel of Fig. [I] do not appreciably contribute 
to the band-gap variation. In fact, this also applies for 
the overall band structure variation. 

Further charge rearrangements 5 that are inher- 
ent to vd W for ces (and require a SC vdW-DF 
calculations) 5 * 52 * 5 -^ are expected to be even smaller than 
those resulting at the the GGA level. The inherent rear- 
rangements are not expected to be of importance for the 
band structure. This justifies our use of non-SC vdW-DF 
in this study of bilayer graphane. 



C. Hybridization and concerted effects on kinetic 
and potential energy 

In the remainder of this section we investigate the role 
of hybridization for the observed band-structure modifi- 
cations. In the simplest picture of hybridization a bond- 
ing and antibonding hybrid WF can be formed when two 



degenerate atomic or layer WFs <pi and <p2 approach each 
other, 



T2 K 



(5) 



If the actual BL WF equals ip+ or alternatively |^bl| 2 = 
|?/; + | 2 , its energy is shifted to lower energies. Similarily, 
if the actual BL WF equals ip~ or alternatively |^bl| 2 = 
|^~| 2 , its energy is shifted to lower energies. The energy 
shifts in such a simple picture of hybridization effects are 
due to a gain and a loss of kinetic energy. 

This simple hydrogen-like picture of hybridization 
needs to be modified in the present case for two reasons. 
First, the unhybridized ML WFs are already complex 
objects possessing internal nodes. Second, the WFs live 
in a background effective potential V. The hybrid WFs 
will then interact with this potential leading to further 
modifications of the actual bonding (B) and actual anti- 
bonding (A) WFs ip B and ip A . The actual hybrid WFs 
ip B and ip A will therefore no longer coincide with ip + and 
ip~ , nor will their energy shift be only of kinetic nature. 

We now move the discussion to a comparison of ML 
and BL graphane. For the ML, we denote the WFs by 
0n,k> where n is the band index and k the wave vector. 
For the BL, ip B k = ^2n-i,k is the bonding WF associated 
with two </> n? k located on different sheets; ^ k = ^2n,k 
is the antibonding WF. 

These WFs (here collectively denoted by </? n ,k) satisfy 
the Kohn-Sham equation 

[-V 2 + (Kff - E F j\ ip nM = (E mM - E F )ip nM . (6) 

Here, V e s is the effective potential (which, in general, is 
different for the ML and BL system), E F is the Fermi 
level (which may also be different in the ML or BL sys- 
tem) and E n ^ the band energy of the WFs (also different 
in general) . Accordingly, we can separate the kinetic- and 
potential-energy shifts of hybrid WFs as 
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In the present analysis we focus on a quantitative eval- 
uation of the per-orbital potential-energy shifts in Eq. (|8| 
and on a qualitative account of the changes in the kinetic- 
energy term in Eq. ([7]). We replace the effective poten- 
tial by the electrostatic potential V^ s (consisting of the 
Hartree potential and the atomic core potentials), ne- 
glecting effects from exchange and correlation, and give 
qualitative accounts of the changes in kinetic energies. A 
quantitative comparison of kinetic-energy shifts would be 
desirable but, since we are using pseudo-WFs, the eval- 
uation of Eq. (JT]) is nontriviaP^and beyond the present 
scope. 

We obtain a qualitative analysis of the changes in ki- 
netic energies by plotting the change of partial electron 
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k 


n 


AE B 


AE A 


AV B 


AV A 








in 


eV 




r 


2 (VB) 


-0.17 


0.14 


0.20 


-0.34 


r 


3 (VB) 


-0.48 


0.51 


0.43 


-1.04 


K 


6 (CB) 


-0.37 


0.31 


0.37 


-0.51 


r 


6 (CB) 


0.00 




-1.07 





TABLE II: Total energies shifts AE and electrostatic con- 
tributions AV according to Eqs. ^ of hybrid WFs in the 
BL with respect to ML WFs at F and K for several bands (n 
specifies the corresponding ML band and CB or VB whether 
this band belongs to the conduction or valence band). 



density associated with an (anti-) bonding BL WF 

p B J£ = \<i A \ 2 0) 

with respect to a sum of or difference between the corre- 
sponding ML WFs, 

P± k H< k | 2 = l/2|<^±0 2 - (10) 
The differences 

PnM ~ PnM ( U ) 

measure the extent to which the BL WFs experience a 
reduction or enhancement of kinetic energy with respect 
to a simple hybridization, Eq. (J5|. 

D. Hybridization in the valence band 




Hybridization in the valence band is found, for exam- 
ple, around the Y point in the bands that correspond to 
band no 2 and 3 in the ML, see Figs. [4] and [5] We em- 
phasize that the hybridization of the corresponding WFs 
is not a signature of binding. The energy splits are (es- 
sentially) symmetric and since both the bonding and the 
antibonding states are occupied, there is no (significant) 
net gain in total energy. 

In Table [nj we list the total-energy shifts AE B / A (ob- 
tained directly from our calculations) and the potential- 
energy contributions to these shifts AV B / A for the cor- 
responding WFs. Interestingly, we find that the electro- 
static contributions to the energy shifts are positive for 
bonding BL WF whereas they are negative for the an- 
tibonding BL WF. The kinetic-energy gain (loss) must 
therefore be significantly larger than the loss (gain) in 
potential energy of the bonding (antibonding) BL WF to 
produce the ordering shown in the bottom panel of Fig. [5] 
(compare states identified by triangles at the T point). 

The top panels of Fig. [6] show sections (through the 
main diagonal of the unit cell) of the difference between 
the partial densities that correspond to the bonding and 
antibonding VB BL WFs at T and the partial densities 
that correspond to the pure sum of (difference between) 
the corresponding ML WFs (with band index 3). The 



FIG. 6: (Color) Qualitative analysis of kinetic-energy con- 
tributions to total-energy shifts of hybrid WFs in terms of 
contour-plot sections of differences between partial electron 
densities of bonding (antibonding) BL WF and su ms o f (dif- 
ferences between) associated ML WFs, Eqs. ^ and |To| . Pos- 
itive (negative) contours are in red (blue). A green color be- 
tween the MLs in the top panel reflects that the BL VB WFs 
at r are essentially simple hybridizations of the ML WFs. In 
the mid-panel, our comparison indicates that the bonding (an- 
tibonding) LCB BL WFs at K possess a lower (higher) kinetic 
energy than the simple bonding-type sum of (antibonding- 
type difference between) ML WFs. The bottom panels show 
that the lowest CB BL WF at T differs significantly from both 
the sum of and the difference between ML WFs and that the 
BL WF possesses a more complex hybridization nature than 
implied by a simple picture. However, the closer resemblance 
to the antibonding difference of ML WFs indicates a loss of 
kinetic energy with respect to an individual ML WF. 



green color (indicating an absence of any decrease or en- 
hancement relative to a simple hybridization) between 
the MLs in both panels shows that the BL WFs gains or 
looses kinetic energy. These kinetic-energy changes evi- 
dently more than make up the concerted changes in the 
potential-energy terms. 
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E. Conduction-band modifications 



For the CB, we focus on the band-gap modifications 
and WFs at the points marked in the top panel of Fig. [5] 
(at K and T). At K, the total-energy splitting leads to the 
reduced band gap. As shown in Tab. [TTJ the potential- 
energy shift is positive for the bonding BL WF at K 
(marked with a downward triangle in Fig.[5| and negative 
for the antibonding BL WF (marked with an upward 
triangle). This observation is in line with those made at 
T in the VB. 

The mid-panels of Fig. [6] show the differences between 
partial density associated with the lowest bonding (next- 
lowest antibonding) CB BL WFs and the density associ- 
ated with the pure sum of (difference between) the cor- 
responding ML WFs at K. For the bonding BL WF, the 
partial density is increased between the two ML with re- 
spect to the pure sum of ML WFs. For the antibonding 
BL WF, the partial density is decreased with respect to 
the difference between the ML WFs. This indicates that 
the kinetic energy gain (loss) of the bonding (antibond- 
ing) BL WF is larger than within a simple hybridization 
picture where the hybrid WFs already possessing a ki- 
netic energy gain (loss) with respect to an individual ML 
WF. 

At T, the contribution of the potential energy to the 
total shift of the lowest lying CB WF (marked with a 
square Fig. [5| is negative, see Table [TlJ In a simple hy- 
bridization picture one would expect a WF with a bond- 
ing nature and an additional decrease of kinetic energy. 
The vanishing shift in total energy, however, requires a 
kinetic-energy offset that compensates for the negative 
potential-energy shift. Specifically, one must therefore 
expect a more complicated hybridization of this BL WF. 

The bottom panels of Fig. [6] show the difference be- 
tween the partial density of the lowest CB BL WF and 
the pure sum of (left panel) and the pure difference be- 
tween (right panel) the corresponding ML WFs at T. The 
significantly negative value of the contours between the 
MLs in the left panel indicates that the BL WF at T has 
a higher kinetic energy than the pure sum of ML WFs. 
Similarly, the kinetic energy is lower than in the differ- 
ence between ML WFs, indicated by the positive value 
of the contours between the MLs in the right panel. 

The contour plots show that the lowest CB BL WF 
at T is not a simple hybridization of ML WFs. Also, 
we notice that the differences in the right panel are not 
as pronounced as the differences in the left panel. This 
suggests that the lowest CB BL WF at T possesses rather 
an antibonding nature (although the lowest CB BL WF 
has a smaller kinetic energy than a simple antibonding 
hybridization). Such an antibonding nature is consistent 
with the vanishing total-energy shift and the negative 
potential-energy shift. 



VI. DISCUSSIONS: BULK GRAPHANE 

The most pronounced electronic difference between 
bulk graphane and ML or BL graphane is the fact that 
the A point takes over the role of the T point, see Fig. [4] 
At A, the band gap is smallest in the bulk and the size of 
that gap coincides essentially with the size of the gap at 
T in the ML or BL. Also, in the bulk, the deviations in 
the band-gap variations can be considerably larger than 
those in the BL, see bottom panel of Fig. [4j Here, we 
give a more detailed analysis of the overall band struc- 
ture effects associated with graphane assembly into bulk. 

In the the left panel of Fig. [7| we show the band dia- 
gram along KT, TM, and MK. Focusing on the CB, the 
bulk band structure is very different from that of the BL 
along the same path, see top panel of Fig. [4] In partic- 
ular at the T point, that gap is approximately twice as 
large as in the BL. Also, the band energy at K is consid- 
erably higher than in the BL. Only at M, we find similar 
band energies. Thus, the effect of vdW-bonding on lo- 
cal features of the band structure and related observable 
properties can be dramatic. 

In the mid-panel of Fig. [7| we show show the band 
diagram along HA, AL, and LH. Again, the the bulk 
band structure is very different from that of the BL along 
the parallel path along KT, TM, and MK. The band 
energy at H is lower than at L, while it is the other way 
around at K and at M in the BL. Only around the A 
point (corresponding to T in the BL and ML cases) do 
the bulk and ML band diagram show similar features. 

In the right panel of Fig. [7| we effectively combine 
the two diagrams from the left and mid-panel, calculat- 
ing the band diagram for a would-be 2-layer unit cell at 
k z = 0. The resulting band diagram shows the features 
of the BL band diagram and of the VB there is also a 
very good quantitative agreement. Therefore, the BL 
band-diagram can partly be understood as a zone-folded 
version of the bulk band-diagram. 



VII. SUMMARY AND CONCLUSION 

This paper predicts and characterizes band-structure 
modifications produced indirectly as a geometry-induced 
effect of pure dispersive binding between two macro- 
molecules in possible new material systems: bilayer and 
bulk grapha ne! 19 I 2Q I Using non-selfconsistent vdW-DF 
calculat ions j 4 * 5 * 1 3 -H^ we determine the binding separa- 
tion (c-lattice parameter) and binding energies in vdW- 
bonded graphane bilayers (bulk). We use the calculated 
separations to obtain corresponding GGA-DFT band di- 
agrams. Our results demonstrate that vdW interactions 
can significantly alter electron behavior, at least locally 
in the Brillouin zone. In graphane, the direct band gap 
is reduced (increased) by up to 1.2 eV (4 eV). 

Our analysis shows that the origin of the band-gap 
modifications in this system is the concerted action of 
two geometry-induced effects. The first cause is the hy- 
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H K,H 



M, L K, H 



FIG. 7: Band diagrams for bulk graphane. The left panel shows the band diagram along the same path in k space as shown 
for the bilayer in the top panel of Fig. [4] The midpanel shows the band diagram along a parallel path in k space with k z = 1/2. 
In the right panel we combine both band diagrams into one plot and reproduce a band diagram that possess essentially all 
features that are present in the BL band diagram, compare top panel of Fig. [4] 



bridization between unoccupied wave functions in the 
lowest conduction band. This effect generally leads 
to lower (higher) energy state for the bonding-type 
(antibonding-type) hybrid wave function, but not always, 
since the graphane wavefunctions can possess a more 
complex hybridization nature than what applies in the 
hydrogen-like case. The second cause is the modified 
electrostatic interaction between the hybrid wave func- 
tions and the electron density. This cause can either 
increase or decrease the energy gain (loss) in a bonding- 
type (antiboning-type) hybridization. Moreover, it may 
result into a modified hybrid wave-function where a pure 
bonding- type (antibonding-type) character is lost. 

Our analysis also shows that selfconsistent charge re- 
arrangements (on the GGA level) with respect to the 
ML density have no significant impact on the band- 
structure. Additional charge rearrangements described 
by self-consistent vdW-DF calculations 5 are expected to 
be even smaller, justifying our use of non-selfconsistent 
vdW-DF. 

The nature of conduction and optical absorption in BL 
graphane would be determined by the region around the 
T point. There, the band structure of the bilayer essen- 
tially coincides with that of the monolayer. Therefore, 
we expect graphane multilayers to behave electronically 
similar to a graphane monolayer, at least for properties 
defined by a simple response. For bulk graphane, the 
behavior is more complicated. We emphasize that other 
vdW-bonded systems may exist where significant band- 
gap modifications arise at Brillouin-zone points having 
higher relevance for the electronic nature of the material. 

Our results for BL graphane suggest that vdW forces 
can have non-negligible indirect effects on the overall 
band structure in layered or macromolecular materials. 
A similar effect can be found in V205p^ where tra- 
ditional GGA severely overestimate the c lattice con- 
stant and where vdW-DF provides a more accurate 



description.^ Furthermore, we notice that vdW bind- 
ing of the intrinsic semimetal graphene to metaP or 
semiconductoiESnj] surfaces seems to generally lead to a 
shift of the Fermi level, rendering graphene a true metal. 
These observations together with the fact that the vdW 
binding strength in surface/adatom systems (and thus 
presumably the strength of the corresponding effect on 
the band structure) depends on the choice of the sub- 
strate material 9 and on the substrate morphology^ im- 
plies a possibility to exploit dispersive interactions also 
in band-gap engineering. 

Finally, a comment on the accuracy of the predicted 
band-gap modifications is in order. GGA DFT typically 
severely underestimates the band gap in semiconductors 
or insulators. A possible remedy is the use of the GW 
method where the Kohn-Sham orbitals are used to con- 
struct the self-energy operator. The GW method also 
contains non-local correlations and therefore provides an 
alternative description of vdW interaction. GW calcula- 
tions come at a considerably higher computational cost 
than that of vdW-DF, however. The procedure illus- 
trated here, using vdW-DF calculations to determine 
binding morphologies followed by band-structure calcula- 
tions might well also be adapted for GW, pursuing char- 
acterizations of band-structure effects in sparse matter. 
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